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ABSTRACT 

We compute from lattice QCD in the valence (quenched) approxima- 
tion the partial decay widths of the lightest scalar glueball to pairs of 
pseudoscalar quark- ant iquark states. These predictions and values ob- 
tained earlier for the scalar glueball's mass are in good agreement with 
the observed properties of /j(1710) and inconsistent with all other ob- 
served meson resonances. 
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It is generally believed that QCD predicts the existence of glueballs, reso- 
nances composed mainly of chromoelectric field without a valence quark-antiquark 
pair, occurring either as physical particles by themselves or in linear combination 
with states which do include a valence quark and antiquark. Whether such states 
have been identified so far in experiment remains ambiguous. A crucial problem is 
that the properties of glueballs are not expected to be drastically different from the 
properties of flavor singlet bosons including valence quarks and antiquarks. Thus 
the identification in experiment of states with large glueball contributions is difficult 
if not impossible in the absence of a reliable evaluation of the properties predicted 
for glueballs by QCD. We believe the lattice formulation of QCD provides the most 
reliable method now available for determining QCD's predictions for the masses and 
decay couplings of hadrons. 

Some time ago we reported a value of 1740(71) MeV for the valence (quenched) 
approximation to the infinite volume continuum limit of lattice QCD predictions for 
the mass of the lightest scalar glueball. This result was obtained using ensembles of 
25000 to 30000 gauge configurations on each of several different lattices. An earlier 
independent valence approximation calculation when extrapolated to the contin- 
uum limit 1^ following Ref. [0], yields 1625(94) MeV for the lightest scalar glueball 
mass. This calculation used several different lattices with ensembles of between 1000 
and 3000 configurations each. If the two mass evaluations are combined, taking into 
account the correlations between their statistical uncertainties arising from a common 
procedure for converting lattice quantities into physical units, the result is 1707(64) 
MeV for the scalar glueball mass. Both the mass prediction with larger statistical 
weight and the combined mass prediction are in good agreement with the mass of 
/j(1710) and are strongly inconsistent with all but /o(1500) among the established 
flavor singlet scalar resonances. For /o(1500) the disagreement is still by more than 
three standard deviations. 

The valence approximation, used in the mass calculation of Refs. [0 , may be 
viewed as replacing the momentum and frequency dependent color dielectric constant 
arising from quark-antiquark vacuum polarization with its zero-momentum, zero- 
frequency limit 1^. This approximation is expected to be fairly reliable for long- 
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distance properties of hadrons. For example, the infinite volume continuum limits 
of the valence approximation to the masses of eight low-lying hadrons composed of 
quarks and antiquarks differ from experiment by amounts ranging up to 6% A 
6% error in the glueball mass would be 100 MeV and, according to an adaptation 
of an argument giving a negative sign for the valence approximation error in f^^ 0, 
the sign of this error is also expected to be negative. Thus the scalar glueball in full 
QCD should lie above the valence approximation mass, and correcting the error in the 
valence approximation should not drastically change the comparison with experiment. 

The most likely interpretation of /o(1500), we beheve, is not as a glueball 
but as a state composed largely of an ss quark-antiquark pair. The su scalar and 
tensor are nearly degenerate at about 1430 MeV. Thus the ss scalar and tensor 
should lie close to each other somewhere above 1430 MeV. Since the ss tensor has 
been identified at 1525 MeV, an ss scalar at 1500 MeV would be quite natural. 

The crucial question not answered by the mass results, however, is whether the 
decay width of the lightest scalar glueball is small enough for this particle actually to 
be identified in experiment. In addition, it is sometimes argued that since glueballs 
are flavor singlets they should have the same couplings to 271q, to 2Kl, and to 2//. 
This expectation is violated by /j(1710) decay couplings. 

In the present article we report the first lattice QCD calculation of the valence 
(quenched) approximation to the partial decay widths of the lightest scalar glueball 
to pairs of pseudoscalar quark-antiquark states. The calculation is done with 10500 
gauge configurations on a single lattice, 16^ x 24, at /3 = 5.70 corresponding to inverse 
lattice spacing = 1.35 GeV. We believe this lattice has spacing sufficiently small 
and volume sufficiently large to give partial widths within 30% of their infinite volume 
continuum limits. The predicted decay couplings, combined with the mass prediction 
of 1740(71) MeV, give a total two-pseudoscalar decay width of 108(29) MeV for the 
scalar glueball. With any reasonable guess concerning the scalar glueball's branching 
fraction to multibody decay modes, the resulting total decay width is well below 200 
MeV and therefore small enough for the scalar glueball to be identified in experiment. 
In fact, the predicted total two-pseudoscalar decay width, and individual couplings 
to 2770, to 2Kl, and to 2ri are all in good agreement with properties of /j(1710) and 
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Figure 1: Decay couplings. 



inconsistent with all other established flavor singlet scalar resonances. A comparison 
of our results with data for /j(1710) is shown in Figure |l]. 

Glueballs found in the valence approximation, according to one simple inter- 
pretation, contain no admixture of configurations with valence quarks or antiquarks. 
Thus we consider the agreement between the mass and decay couplings found in the 
valence approximation and the observed mass and decay couplings of /j(1710) to be 
strong evidence that this state is largely a scalar glueball with at most some relatively 
smaller amplitude for configurations including valence quark-antiquark pairs. 

The calculations presented here were carried out on the GFll parallel computer 
[y] at IBM Research and took approximately two years to complete at a sustained 
computation rate of between 6 and 7 Gfiops. A preliminary version of this work is 
discussed in Ref. [|10 |. 

In the remainder of this paper we describe our method for determining scalar 
glueball decay couplings then present our numerical results. 

To evaluate glueball decay couplings we work with a euclidean lattice gauge 
theory, on a lattice x T, with the plaquette action for the gauge field, and the 
Wilson action for quarks. It is convenient initially to assume exact flavor SU(3) 
symmetry for the quark mass matrix. With each gauge configuration fixed to lattice 
Coulomb gauge, we construct a collection of smeared fields. We describe smearing 
only for the particular choice of parameters actually used in the decay evaluation. 
Let Ui{x) for a space direction i = 1,2,3, be a smeared link field ^ given by the 
average of the 9 links in direction i from the sites of the (3 site) x (3 site) square 
oriented in the two positive space directions orthogonal to i starting at site x. Let 
Vij{x) be the trace of the product around the outside of a (3 link) x (3 link) square 
tr[Ui{x)Uj{x + 3i + 2j)Uj{x + 5j)U}{x-2i + 2j)], where i is an z-direction unit vector. 
Define the zero-momentum scalar glueball operator g{t) to be the sum of the Vij{x) 
for all i,j and x with time component t. Let the quark and antiquark fields \E'(x) 
and "^{x) be Wilson quark and antiquark fields smeared by convoluting the local 
Wilson fields with a space direction gaussian, invariant under lattice rotations and 
with mean-square radius 6.0. The smeared pseudoscalar field 7ii{x) with flavor index 
i is \E'(x)7^Aj\E'(a;), where Aj is a Gell-Mann flavor matrix. Let 7fj(A;,t) be the Fourier 
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transform of ni^x) on the time t lattice hyperplane. 

Define and i^J be the energy of a single pseudoscalar at rest or with 
momentum magnitude \k\ = 27r/L, respectively. The field strength renormalization 
constant rj'^ is defined by the requirement that for large t the vacuum expectation 
vahie < 7f,|(0, t)7r,;(0, 0) > approaches {rijyL'^exp[—Eit]. Define similarly from a 
pseudoscalar field with momentum magnitude = 27r/L. In the valence approxi- 
mation, the glueball is stable so that its mass and field strength renormalization 
constant rj^ can be defined by the requirement that, for large t, < g{t)g{0)) > ap- 
proaches {r]^)'^L^exp{—E^t). 

Prom pseudoscalar fields at position and times ti, define the two-pseudoscalar, 
flavor singlet fleld 11(^1,^2) to be (16)~^/^ Xii 7ri(0, ti)7ri(0, ^2), where the sum over 
i runs from 1 to 8. Let the zero-momentum, two-pseudoscalar flavor singlet fleld 
ni(ti,i2) be (16)"^/^ 7fi(0, ti)'n-i(0, ^2)- Deflne the two-pseudoscalar fleld 62(^1,^2) 
to be (24)~^/^ X^.g7fi(^, ti)7fi(— ^, t2) where the sum for k is over the three positive 
orientations with |^| = 27:/ L. 

Let |1 > and |2 > be, respectively, the lowest and second lowest energy 
flavor singlet, rotationally invariant two-pseudoscalar states. Both states are nor- 
mahzed to 1. Let Ef^ be the energy of \i >. Deflne the amphtudes r]fj^{t) to be 
L-^ < i\Uj{t,0)\n >. For large t, rj^fit) has the asymptotic form r]^^exp{-E]t). 
The diagonal coefficients rfl^ and TyfJ are expected to be larger than the off-diagonal 
ryjf and 77]^^, respectively. As a consequence of the interaction between pairs of pseu- 
doscalars, however, the off-diagonal coefficients will not be zero. 

Connected three-point functions from which coupling constants can be ex- 
tracted are now given by Ti(tc,,t^) deffned as < g{tg)Ili{t.,^,Q) > — < g{tg) >< 
^iitn, 0) >. If the quark mass, and thus the pseudoscalar mass, is chosen so that E^'^ 
is equal to E^, the lightest intermediate state which can appear between the glueball 
and pseudoscalars in a transfer matrix expression for Ti(t^,t^) is |1 >. Thus for large 
enough tg with fixed, Ti(tg,t^) will be proportional to the coupling constant of a 
glueball to two pseudoscalars at rest. If the quark mass is chosen so that i^J'^ is equal 
to E^, however, the lightest intermediate state which can appear between the glueball 
and pseudoscalars in a transfer matrix expression for T2{tg,tT^) we still expect to be 
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|1 >, not |2 >, since is expected not to be zero. To obtain from T2(tg,t,r) the 

coupling of a glueball to two pseudoscalars with momenta of magnitude 27rL~^, the 
contribution to T2{tg,tT^) arising from the |1 > intermediate state must be removed. 
From the three-point functions we therefore define the amphtudes 

Si{tg,t-^) — Ti(tg,t-^) — -—j—^Tj{tg,tT^), (1) 

for {i-ij) of either (1,2) or (2,1). In S2{tg,tT^) the contribution of the undesirable 
|1 > intermediate state has been canceled. In Si{tg,tT^) a contribution from the 
intermediate state |2 > has been canceled. Although the subtraction in Si{tg,tT,) is 
irrelevant for large enough tg, we expect that as a result of this subtraction Si{tg, i^) 
will approach its large tg behavior more rapidly than does Ti(ig,i^). 

An additional intermediate state which can also appear in a transfer matrix 
expression for either Tj(tc,,t^) is the isosinglet scalar bound state of a quark and an 
antiquark. For the parameter values used in the present calculation we have found 
that this state has a mass in lattice units above 1.25 while the scalar glueball mass is 
0.972(44). Thus for large enough tg the scalar quark-antiquark state will make only 
its appropriate virtual contribution and does not require an additional correction. 

At large tg and t-^., the three-point functions become 

Si{tg,t^) ^ / = Si{tg,tTr), (2) 

where ci = l/\/2, C2 = \/3, r is (»7r2''^2f )/(^ri'^22') Ai and A2 are the glueball 
coupling constants to a pair of pseudoscalars at rest or with momenta of magnitude 
27: L~^, respectively. The factors rjfj^ are given by the large t behavior of r]fj^{t) as 
discussed earlier. For T ^ tg > tj^, the factors Si{tg,tT^) are 

Si{tg,t^) = J2^M-E'\t-tg\-E^\t\- 

t 

E^\t-t^\-5iit,t^)\t-t^\], (3) 

where, for t > ^tt, 5i(^,^7r) is the binding energy E^'" — 2EJ and otherwise Siit^tT^) is 
0. 
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The coupling constants in Eq. (|^) have been identified by comparing Si(tg, t^,) 
with the three-point functions arising from a simple phenomenological interaction 
lagrangian. This procedure is correct to leading order in the coupling constants. A 
similar relation used to find coupling constants among hadrons containing quarks has 
recently yielded several predictions in good agreement with experiment |]ll|. The Aj 



are normalized so that in the continuum limit they become, up to a factor of — i, 
Lorentz-invariant decay amplitudes with the standard normalization convention used 
in the section on kinematics of the Review of Particle Properties. 

To obtain values of Aj from Eq. @) we need the amplitudes rj^-'it). These 
we determine from propagators for two-pseudoscalar states. Define two-pseudoscalar 
propagators Cj(ti,t2) to be < Ii{ti + 2^2,^1 + ^2)1^4(^2,0) >. For moderately large 
values of ti, these amplitudes approach 

Ci{hM) = Cuexp{-Erh)+C2^exp{-Erh), (4) 

a, = Vait2Krit2) + V6v^''{t2Krit2). (5) 

From these expressions the required rifj^{t) can be extracted. 

The rjfj'it) in Eq. @ serve, among other purposes, to correct for the interaction 
between the two pseudoscalars produced by a glueball decay. In the valence approx- 
imation this interaction does not include the production and annihilation of virtual 
quark- antiquark pairs. Correspondingly, in the numerical evaluation of 6*^(^1,^2) from 
quark propagators, we include only terms in which all initial quarks and antiquarks 
propagate through to some final quark or antiquark. Terms in the two-pseudoscalar 
propagator in which initial quarks propagate to initial antiquarks can be shown to 
correspond to processes missing from glueball decay in the valence approximation. 
For very large ti and T, the Cj(ti, ^2) are given by a sum of two terms each of which 
is a slightly more complicated version of one of the exponentials in Eq. (^). This 
complication occurs, for example, because in the valence approximation the exchange 
of a p between the pseudoscalars produced in a glueball decay is not iterated in 
the same way as in full QCD. Each term in Eq. @ holds without modification if 
I^TTTT _ 2^;^'^ 12^2^2 << 1. The intervals of ti we use to determine the r]ff fall well 
within this limitation. In any case, as we will discuss below, the measured values of 
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7]'^-' turn out to be close to their values for noninteracting pseudoscalars. As a con- 
sequence, the corrections due to interations between the decay pseudoscalars which 
the f]'^-' contribute to the predicted values of Aj are comparatively small. 

We now turn to our numerical results. At /3 = 5.7 on a 16^ x 24 lattice, with 
an ensemble of 10500 independent configurations, we determined glueball and single 
pseudoscalar energies and renormalization constants following Refs. and re- 
spectively. For , as mentioned above, we found 0.972 ± 0.044. On a lattice of size 
16^ X 40 we then evaluated the two-pseudoscalar propagator Cj(ti,t2) at k = 0.1650 
using 100 independent configuration, and at k = 0.1675 using 875 independent con- 
figurations. Fitting the ti dependence of Cj(ti,t2) to Eqs. (Q) and (|^), we determined 
-Ef^ and Vii^{t2) for a range of different t2- At k = 0.1650 we obtained results for 

< t2 < 4, and at 0.1675 we found results for < t2 < 5. The values of E^'^ 
were statistically consistent with being independent of t2 in all cases. The i]^j^{t2) 
were consistent with the asymptotic form ri^J^exp{—EJt2) in all cases for ^2 ^ 2. At 
K = 0.1650 for E^"^ we obtained 0.908(5), giving glueball decay to |1 > nearly on 
mass shell. At K = 0.1675 for E^^ we found O.SdStomi 

giving glueball decay to 

|2 > nearly on mass shell. For the normalized ratios fj^j^ defined as Tjfj^ /{rij)'^, at 
K = 0.1650 we obtained for ij of 11, 12, 21 and 22, the values 0.988(30), 0.091(8), 
-0.087(8), and 1.065(13), respectively. At k = 0.1675 we found 1.050(21), 0.107(6), 
-0.112(8), 1.053(53). For noninteracting pseudoscalars fjfj^ is 1 for i = j and other- 
wise. Our data is close to these values. The final value of Ai is changed by less than 

1 standard deviation and the final A2 is changed by less than 2 standard deviations 
if we ignore the determination of fj^J^ and simply use the the noninteracting values. 

From our 10500 configuration ensemble on a 16^ x 24 lattice, we evaluated 5*1 
and 5*2 for glueball decay on mass shell at k of 0.1650 and 0.1675, respectively. We 
obtained statistically significant results for < tg — t,, < 2 with < < S. At each 
point within this range we then determined effective Aj using Eq. (^). We found Ai 
and A2 statistically consistent with being constant for t,, > 3 and t,, > 2, respectively, 
and all values of tg—tT^. Figure H, for example, shows effective A2 in units of the p mass 
as a function of for tg — t^^ = 2, in comparison to a fit with 2 < t^^ < Q, tg — = 2. 
Figure ^ shows fitted values of A2 on the interval 2 < < 6 for fixed tg — t^, of 0, 
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Figure 2: X2 ioi tg — t^^ — 2. 
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Figure 3: A2 fitted on 2 < < 6 as a function of tg — t^^ for the fitting interval. 
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1 or 2. To extract final values of Aj, we tried fits to all rectangular intervals of data 
including at least 4 values of and at least 2 values of tg — t-j^- For each Aj we chose 
the fit giving the lowest value of per degree of freedom. The window determined 
in this way for Ai is 3 < < 7 with 1 < — tjr < 2, and for A2 is 2 < < 6 with 
< tg — t^^ < 1. The horizontal line in Figure ^ shows the final value of A2. Over 
the full collection of windows we examined, the fitted results varied from our final 
results by at most 1 standard deviation. We believe our best fits provide reasonable 
estimates of the asymptotic coefficients in Eq. (0). 

So far our discussion has been restricted to QCD with u, d and s quark masses 
degenerate. An expansion to first order in the quark mass matrix taken around some 
relatively heavy SU(3) symmetric point gives glueball decay couplings for tt's, K's and 
r^'s which are a common linear function of each meson's average quark mass. Since 
meson masses squared are also nearly a linear function of average quark mass, the 
decay couplings are a linear function of meson masses squared. Thus from a linear 
fit to our predictions for decay couplings as a function of pseudoscalar mass squared 
at unphysical degenerate values of quark masses we can extrapolate decay couplings 
for physical nondegenerate values of quark masses. From this linear fit a prediction 
can also be made for the decay coupling of the scalar glueball to 77 + r]', if we ignore 
the contribution to the decay from the process in which the 77 quark and antiquark 
are connected to each other by one propagator and the 1]' quark and antiquark are 
connected to each other by a second propagator. 

Figure |I] shows predicted coupling constants as a fuction of predicted meson 
mass squared along with linear extrapolations of the predicted values to the physical 
TT, K and rj masses, in comparison to observed decay couplings^] for decays of /j(1710) 
to pairs of vr's, K's and 77's. Masses and decay constants are shown in units of the p 
mass. Our predicted width for the scalar glueball decay to r] + rj' is 6(3) MeV. For 
the ratio A^^//A^^ we get 0.52(13). We predict a total width for glueball decay to 
pseudoscalar pairs of 108(28) MeV, in comparison to 99(15) MeV for /j(1710). 
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